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Abstract. We study the chemodynamical evolution of elliptical galaxies and their 
X-ray and optical properties using high-resolution cosmological simulations. Our 
Tree N-body/SPH code includes a self-consistent treatment of radiative cooling, 
star formation, supernovae feedback, and chemical enrichment. We present a series 
of ACDM cosmological simulations which trace the spatial and temporal evolution 
of abundances of heavy elements in both the stellar and gas components of galax- 
ies. A giant elliptical galaxy formed in one of the simulations is quantitatively 
compared with the observational data in both the X-ray and optical regime. We 
implement a treatment of both pre-heating and AGN heating in this simulation, 
and examine the effect of these processes on elliptical galaxy formation. We find 
that the adopted pre-heating (T = 10^ K at 2; = 4) is not strong enough to explain 
the observed X-ray or optical properties. On the other hand, our AGN heating 
model in which the gas inflow induces the AGN heating is consistent with both 
the X-ray and optical properties. 
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1. Introduction 

Developments in both ground- and space- 
based observational facilities has made it 
possible for astronomers to understand de- 
tailed self-consistent optical. X-ray, and 
theoretical analyses of galaxy formation 
and evolution. For the case of elliptical 
galaxies, optical observations provide con- 
straints on the properties of the stellar com- 
ponent, while X-ray observations constrain 
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the physical conditions of the hot interstel- 
lar medium (ISM). 

iKawata CibsTTTI EfMh (KB03b) 
presented a first attempt to explain both 
the X-ray and optical properties of ob- 
served elliptical galaxies via the use of self- 
consistent cosmological simulations. Using 
a standard "recipe" for galaxy formation, 
they found that radiative cooling is impor- 
tant to interpret the observed X-ray lumi- 
nosity (Lx), temperature (Tx), and metal- 
licity ([Fe/H]x) of the hot gas of elliptical 
galaxies. However, they also found a seri- 
ous problem in that the cooled gas also 
leads to excessive star formation at low 
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redshift, and therefore results in underly- 
ing galactic stellar populations which are 
too blue with respect to observations. Their 
study demonstrates that the cross-check 
over both X-ray and optical properties pro- 
vides stronger constraints for the theoreti- 
cal models (more so than anyone taken in 
isolation), and is essential for any successful 
scenario of elliptical galaxy formation and 
evolution. 

KG03b suggested that a heating mech- 
anism which was not included in their 
simulations would be required to suppress 
this enhanced cooling and consequent star 
formation. In this paper, we examine two 
such potential mechanisms: pre-heating 
and heating by an Active Galactic Nucleus 
(AGN). Previous analytical and numerical 
studies have shown that the prc-heating of 
gas at high redshift may help to explain 
the X-ray luminosi ty and temperature o f 
the hot gas (e.g. iBorgani et al. I l2002() . 
although the source of the pre-heating 
remains unclear. On the other hand, recent 
observations suggest that all elliptical 
galaxies have a cen t ral b lack hole (e.g. 
iFerrarese fc Merritt I 120001) . One might 
expect therefore that once the cold gas 
falls into the central part of a galaxy, the 
AGN may become "active" and heat up 
the surrounding gas, and might be able to 
blow out the gas from the system via a 
radio jet. Once the gas accretion is halted 
by such AGN heating, the AGN itself 
loses its fuel "source" and goes back to a 
quiescent state, which allows the gas to 
cool again. This self-regulation cycle may 
be capable of suppressing star formation 
for a significant period of time, without 
changing the effect of radiative cooling, 
which is in turn required to consider the 
indication for the observed X-ray proper- 
ties. Several earlier studies of the effect of 
AGN heating on the hot gas of elliptical 
galaxi es can be seen in [B inncy & Tabor ^ 
1119951): iBrighenti fc M^hews I |20Q^;' 
iBriiggen & Kaiserl ^^2002^ . Our current 
study examines the effect of pre-heating 
and AGN heating on the optical properties 
as well as the X-ray properties in the 



context of a more self-consistent treatment 
of cosmological evolution. 

2. Methods 

Our simulations were carried out using 
GCD+, our original galactic chemodynam- 
ical evolution code . Details of the code 
are presented in Kawatal l)l999D and 
iKawata fc Gibson I l|2003albD . In GCD+, the 
dynamics of coUisionless dark matter and 
stars is calculated using a gravitational 
Tree N-body code, and the gas compo- 
nent is modeled using Smoothed Particle 
Hydrodynamics (SPH). We calculate ra- 
diative cooling, star formation, chemical 
enrichment, and supernovae (SNe) feed- 
back, self-consistently, and take into ac- 
count both Type la and Type II SNe. We 
assume that SNe feedback is released as 
thermal energy. 

We have carried out a series of 
high-resolution simulations within a stan- 
dard ACDM cosmology (f7o=0.3, Ao=0.7, 
Qh^Omdh-^, h=0.7, and cr8=0.9). Gas 
dynamics and star formation are included 
only within the relevant high-resolution re- 
gion ('^12 Mpc at z=Q); the surround- 
ing low-resolution region ('^43 Mpc) con- 
tributes to the high-resolution region only 
through gravity. The mass of individual gas 
particles in the high-resolution region was 
5.9 X 10'' M0. We identified an appropri- 
ate elliptical galaxy analog in the high- 
resolution region, which acts as the focus 
for this study. The total virial mass of this 
target galaxy is 2 x lO^'^ M©, similar in 
size to that of NGC 4472, a bright ellip- 
tical galaxy in the Virgo Cluster. The tar- 
get galaxy is relatively isolated, with only 
a few low-mass satellites remaining at z=0. 
This target galaxy is the same system dis- 
cussed in KG03b. Fig. 1 of KG03b shows 
the morphological evolution of dark mat- 
ter in the simulation volume, and the evo- 
lution of the stellar component of the tar- 
get galaxy. The galaxy forms through con- 
ventional hierarchical clustering between 
redshifts z=3 and z=l; the morphology 
has not changed dramatically since z=l. 
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Fig. 1. Comparison of the simulated and observed (crosses with error bars) 
Lx — Tx relations {left panel) and [Fe/H]x — Lx relations {right panel). The cir- 
cle/square/diamond/triangle indicates the predictions of Model 1/2/3/4. Open stars in 
the right panel show the predicted mean metallicity of the stellar component for th e four 
models. The observational data are from lMatsushita. Ohashi. fc Makishima I l)200nl) . The 
solid star presents the position of NGC 4472. 



KG03b presented the results of three dif- 
ferent radiative cooling and SNe feedback 
models. We use two of their models as ref- 
erence models for this study; Model 1 here 
corresponds to Model B of KGOBb, which 
includes cooling and weak feedback; Model 
2 here corresponds to Model C of KG03b 
which mimics Model 1, but incorporates 
100 times greater thermal energy (10^^ erg) 
per supernova than Model 1. Since in this 
study we are interested in the effect of the 
pre-heating and the heating by AGN, we 
construct the additional two models: 

Model 3: pre-heating model. The gas 
component in the entire simulation region 
is heated to a temperature of T = 10^ K at 
redshift z = 4. Weak thermal SNe feedback 
(as in Model 1) is adopted. 

Model 4: AGN heating model. The 
most bounded star particle in the target 
galaxy at z = 1 is assumed to be the heat- 
ing source, i.e. the AGN "particle". Since 
the AGN heating requires gas fueling, we 
assume that the AGN heating is active 
only when the divergence of the velocity 
field of the surrounding gas of the AGN 
particle is negative, (V • t>) < 0. The diver- 



gence of the velocity field of the neighbour 
gas particles is calculated using the SPH 
scheme. A constant thermal energy of 
10^'^erg s~^ is deposited to the neighbour 
gas particles when the above condition 
is satisfied. This energy is roughly con- 
sistent with observat ional estim ate s (e.g . 
iBlanton. Sarazin fc M cNamaxal l2003t) . 
Strong thermal SNe feedback (as in Model 
2) is also adopted. 

For all the models, we examine both 
the resulting X-ray and optical properties, 
comparing them quantitatively with obser- 
vation (see KG03b for details). The gas par- 
ticles in our simulations carry with them 
knowledge of the density, temperature, 
and abundances of various heavy elements. 
Using the XSPEC vmekal plasma model, 
we derive the X-ray spectrum for each gas 
particle, and synthesize them within the 
assumed aperture (R'^35 kpc). We next 
generate "fake" spectra with the response 
function of the XMM EPN detector, as- 
suming an exposure time (40 ks) and tar- 
get galaxy distance (17 Mpc). Finally, our 
XPSEC fitting provides the X-ray weighted 
temperatures and abundances of various el- 
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Fig. 2. Comparison of the simu- 
lated f7 - y and y - X CMRs (cir- 
cle/square/diamond/triangle for Model 
1/2/3/4) and those of the Coma cluster 
ellipticals (c rosses) . The observatio n al dat a 
are from iBower. Lucev. fc Ellis I l|l992() . 
The dashed line shows the OMR fitted to 
the Coma Cluster galaxies. The solid star 
shows the position of NGC 4472. 



ements. Conversely, the simulated star par- 
ticles each carry their own age and metal- 
hcity "tag", which enables us to generate 
an optical-to-near infrared spectral energy 
distribution for the target galaxy, when 
combined with our population synthesis 
code adopting simple stellar populations of 
iKodama fc Arimoto 1 1)1997|) . 



3. Results and Conclusion 

To study both X-ray and optical properties, 
we examine the Lx — Tx and Lx — [Fe/H]x 




t (Gyr) 



Fig. 3. Time variation of the star forma- 
tion rate for Models 2 (upper), 3 (middle), 
and 4 (lower). 



relations for the X-ray properties, and the 
colour-magnitude relation (CMR) for the 
optical properties. Figs. ^ and El show the 
Lx — Tx and Lx — [Fe/H]x relations and 
the CMRs for the four models, and com- 
pare the model results with the observa- 
tional data. Models 1 and 2 are re-plot 
from Models B and C of KG03b. As shown 
in KG03b, current standard galaxy for- 
mation models, such as Models 1 and 2, 
can roughly reproduce the Lx — Tx and 
Lx — [Fe/H]x relations. However, the op- 
tical colours of the resulting stellar compo- 
nents for these models are inconsistent with 
the observational data (Fig. ^ , due to ex- 
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Fig. 4. Density vs radius (left) and density 
vs temperature (right) distributions of gas 
particles for Models 1 (upper), 2 (middle), 
and 3 (lower). The solid (dotted) curves 
separate the region where the cooling time 
is shorter (upper region) and longer than 
the Hubble (dynamical) time. 



cessive star formation at low redshift (Fig. 
EJ, which is induced by radiative cooling. 

The pre-heating model. Model 3, leads 
to too high X-ray luminosity and too low 
temperature (Fig. P). Fig. 0] shows that 
in Model 2 the gas whose cooling time, 
tcooi, is shorter than the Hubble time, to, 
is changed into the cold gas. However, 
in Model 3 the gas whose tcooi is longer 
than the local dynamical time, tdyn — 
{Sn/lGGpy^^ , can stay in the hot phase, 
which leads to the higher hot gas density, 
and thus the higher Lx- Also Fig. |2| shows 
that the pre-heating cannot suppress star 
formation until z = 0, although star for- 
mat ion is suppressed unti l z ~ 0.5 (see 
also iTornatore et al. Il2003|) . Consequently, 
a larger amount of remaining gas leads to a 
higher star formation rate at z = 0, which 
results in a bluer colour than Models 1 



and 2. This result demonstrates that the 
pre-heating we assumed here is not strong 
enough for this elliptical galaxy. A higher 
amount of the pre-heating will be exa mined 
in a future work (see also .Boreani et al. I 
l2n?)l . 

The most successful model is the AGN 
heating model, Model 4. Although the X- 
ray luminosity is slightly lower compared 
with the observed value of NGC 4472 (Fig. 
^1, Model 4 roughly reproduces both the 
X-ray and optical observational data. Fig. 
4 shows that our AGN heating model gets 
rid of the cold gas in the central region, but 
does not change the hot gas profile signif- 
icantly. As a result, the hot gas properties 
for Model 4 are similar to Model 2. This 
is due to the self-regulated AGN heating 
which guarantees that the hot gas is not 
"over-heated". We also found that, com- 
pared with Model 2, a larger fraction of iron 
ejected from stars is blown out from the 
system by the AGN heating, which results 
in a lower metallicity in the hot gas (Fig. 
^1. In addition, this AGN heating is pow- 
erful enough to stop the recent star forma- 
tion, and leads to the observed red colours 
of elliptical galaxies (Fig. 

Our cosmological chemodynamical code 
makes it possible to undertake quantita- 
tive comparisons between numerical simu- 
lations and observational data in both the 
X-ray and optical regime with minimal as- 
sumptions. We find that the AGN heating 
induced by cold gas inflow can explain both 
the X-ray and optical observations of el- 
liptical galaxies. Although our AGN model 
is rather simple and we focus on only one 
elliptical galaxy model, this result encour- 
ages more serious studies of the effect of 
the AGN heating on elliptical galaxy for- 
mation. 
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